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Abstract. 

Parton propagation in dense nuclear matter results in elastic, inelastic and coherent 
multiple soft scattering with the in-medium color charges. Such scattering leads to calculable 
modifications of the hadron production cross section that is evaluated in the framework of the 
perturbative QCD factorization approach. Final state medium-induced gluon bremsstrahlung is 
arguably the most efficient way of suppressing large transverse momentum particle production in 
nucleus-nucleus collisions. The observed hadronic attenuation, known as jet quenching, can be 
related to the properties of the medium, such as density and temperature, and carries valuable 
information about the early stages of heavy ion reactions. Non-Abelian energy loss in the quark- 
gluon plasma can be studied in much greater detail through the modification of the two particle 
back-to-back correlations. Perturbative calculations give good description of the redistribution 
of the lost energy in lower transverse momentum particles and predict significant increase of the 
correlation width of away-side di-hadrons. In contrast, energy loss in cold nuclear matter was 
found to be small but for large values of Feynman-x is expected to complement the dynamical 
higher twist shadowing in experimentally observable forward rapidity hadron suppression. 


1. Energy loss in hot and dense nuclear matter 

Jet quenching, the suppression of high transverse momentum hadron production relative to 
the expectation from p+p collisions scaled by the number of elementary nucleon-nucleon 
interactions, has been predicted ^ to be one of the key signatures for the creation of a hot 
and dense quark-gluon plasma (QGP) in ultra-relativistic nuclear collisions. Since, the non- 
Abelian energy loss of fast quarks and gluons in dense nuclear matter has been studied in great 
detail in a variety of theoretical approaches |21 ■ Its key features can be illustrated on the example 
of the reaction operator formalism |3j, which systematically expands energy loss in terms of the 
correlations between multiple scattering centers. 

For A-|-A collisions, the dominant final state energy loss follows the hard partonic scattering, 
which takes place on a time scale to ~ 1/Q- Even in the absence of a medium, the large virtuality 
leads to non-Abelian bremsstrahlung. For real gluons of small and moderate frequencies oj [Sj 
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In Eq. 0* = arcsin(A;_L/t(;) is the angle relative to the jet axis, 8 is the azimuthal cone angle 
(both illustrated in right panel of Eig.^), ag is the strong coupling constant and Cr = 3 (4/3) for 
gluon (quark) jets, respectively. Virtual gluon corrections remove the cu —> 0 infrared singularity 
in the cross sections in accord with the Kinoshita-Lee-Nauenberg theorem but the collinear 
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Figure 1. Left panel: diagrams contributing to the lowest order in opacity energy loss 
expansion. Right panel: the angular distribution of medium-induced bremsstrahlung of 
E = G GeV gluon jet for fixed values of the radiative gluon energy a; = 0.5, 1, 2, 4 GeV. Top 
and bottom panels represent (1+1)D Bjorken expanding medium of transverse size L = 9 fm 
and L = 3 fm, respectively. Inserts show {9*) versus to. 


0* —> 0 divergence has to be regulated or subtracted in the parton distribution functions (PDFs) 
and the fragmentation functions (FFs). 

In contrast, the final state medium-induced bremsstrahlung spectrum is both collinear and 
infrared safe. To first order in the mean number of soft interactions in the plasma the distribution 
of gluons in angle and frequency can be written as 
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In Eq. m a = Z(/c±,g_L), ^g{z) is the position-dependent gluon mean free path and L is 
the transverse size of the medium. The momentum transfers q± are distributed according 
to the normalized elastic scattering cross section (T~i^daei/d?q±{z) = fi‘jj{z)/7r{q‘j_ + iJ,‘jj{z))‘^. 
In this model {q‘j_) oc and for a quark-gluon plasma in local thermal equilibrium 

dvra^T^. For the case of (1-|-1)D dynamical Bjorken expansion of the QGP 
Hd{z) = ^,d{zo) {zo/^g{z) = ^g{zo) {z/zo)^^^- At small angles 9* both the direct and 
double Born terms jH] in Eq. m are divergent. However, they come with different signs 
oc (—= —1, oc i{—i) = 1, equivalent to a phase difference vr. The soft and collinear divergences 
then exactly cancel jHj. 

From Eq. © the gluon distribution is not only hnite when 0* —> 0 but vanishes on average due 
to the uniform angular distribution of momentum transfers from the medium, da cos a = 0. 
We have checked that for physical gluons of k± < uj the cancellations discussed here persist to 
all orders in the mean number of scatterings [Hj. The small frequency and small angle spectral 
behavior of dN^^^/diod sin 9*d6 remains under perturbative control. Given the vastly different 





































angular behavior of the vacuum and the medium-induced gluon bremsstrahlung, Eqs. m and 
it is critical to identify the phase space where cancellation of the color currents induced by 
the hard and soft scattering occurs. We fix the parameters of the medium to ^inizo) = 1-5 GeV 
and Xg{zo) = 0.75 fm at initial formation time zq = 0.25 fm. Since small frequency emission 
is suppressed, we use only a moderate ag = 0.25. Triggering on high pTi hadron directs its 
parent parton “c” away from the medium and places the collision point of the lowest order (LO) 
ab —> cd underlying perturbative process [HI close to the periphery of the nuclear overlap region. 
Then, it is the back-scattered jet “d” that traverses the QGP. For large nuclei, such as Au and 
Pb, path lengths L = 9 (3) fm are used to illustrate central (peripheral) collisions, respectively. 
We limit gluon emission to the forward jet hemisphere, 0 < 9* < ^. 

The angular distribution of medium-induced radiation for F = 6 GeV gluon jet for select 
values of uj is shown in Fig. ^ We find that gluon emission is strongly suppressed within a cone 
of opening angle 6* ~ 0.25 rad. The broad gluon distribution can be characterized by the mean 
emission angle 
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given in the insert of Fig. Q We emphasize that the angular spectrum is distinctly non-Gaussian 
and the mechanism that determines its form is based on the destructive interference of color 
currents from hard and soft parton scattering rather than a slow random walk of the gluon in 9*. 
The large angle distribution or radiative quanta is evident for both large and small frequencies uj. 

In the case of heavy quark creation and propagation the vacuum [Jj and the Bertsch-Gunion 
radiation [S] are modified by the large mass Mq. The mass enters in the propagators as follows, 
—> k‘j_ + Mq /, and suppresses the gluon bremsstrahlung for sin0* < Mq/E. Such 
depletion can only be effective if gluon radiation peaks at 9* —> 0. Since for light quarks 
the final state medium-induced non-Abelian bremsstrahlung is not dominated by small angle 
emission, see left panel of Fig. ^ effects other than the “dead cone” control the energy loss 
pattern of heavy versus light quarks. 

Qualitatively, for light quarks the behavior of the energy loss as a function of the density and 
the size of the system can be summarized to first order in opacity as follows: 
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Eq. 0 assumes energetic jets and neglects kinematic bounds. Here, A_i_ is the transverse size 
of nuclear matter. For static systems AE depends quadratically on the nuclear size. For the 
case of longitudinal Bjorken expansion this dependence is reduced to linear [^j but the energy 
loss is sensitive to the initial parton rapidity density dN^/dy. Understanding the effective color 
charge density dependence of AE the is the key to jet tomography [Hj. 

The calculated energy loss can be incorporated in the perturbative QGD factorization 
approach as a modification of the single and double inclusive hadron production cross sections. 
To lowest order in the double collinear limit these are given by: 
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Figure 2. Left panel: nuclear modification factor for single inclusive pion production 

at SPS, RHIC and the LHC. Right panel: the same nuclear modification for the RHIC Au+Au 
run at c.m. energy of 62 AGeV. 
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In Eq. © iL = 2 is a next-to-leading order A-factor, Xa,b = Pa,b/PNa,Nb the momentum 
fractions of the incoming partons and zi ^2 = Phi,h 2 /Pc,d are the momentum fractions of the 
hadronic fragments. We use standard lowest order Gluck-Reya-Vogt PDFs m and Binnewies- 
Kniehl-Kramer FFs HH. Renormalization, factorization and fragmentation scales are suppressed 
everywhere for clarity. The spin (polarization) and color averaged matrix elements \Mab^cd\'^ 
are given in [^. 

Dynamical nuclear effects in multi-particle production can be studied through the ratio m 
j^{n) ^ da%'^^/dyi • • • dynd^pr^ • • • d'^pr^ 

{^ab) /dy\--- dynd?pTi ■ • • d?pT^ 

Centrality dependence is implicit in Eq. m and the modified cross section per average 
collision da ^^'^"/can be calculated from Eqs. (@J and 0 including the numerically 
evaluated energy loss. Calculation of the nuclear modification in central Au+Au collisions [2] 
at y/SNN = 17,200,5500 GeV are shown in the left panel of Fig. 0. The right panel of 
Fig. 0 extends the GLV jet quenching predictions m to the intermediate RHIC energy of 
y/SNN = 62 GeV. In this case parton energy loss, driven by the medium density dN^/dy, is 
smaller but its effect at high pT is amplified by the steepness of the underlying partonic spectra. 
Similar results for the light pion suppression were found in m. Baryons, however, are expected 
to be significantly less suppressed. In fact, the p/'K~^ ratio was predicted to be larger than at the 
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Figure 3. Left panel: nuclear modification central and peripheral Au+Au 

collisions at RHIC of the back-to-back di-hadron correlations with and without the contribution 
of medium-induced bremsstrahlung. Right panel: angular correlations of hadrons dominated 
by large angle gluon emission. Solid and dashed lines give separately the contribution of the 
radiative gluons and the attenuated parent parton. 


top RHIC energy m Relating the observed high px hadron attenuation ununmais] to the 
energy density in central Au-|-Au collisions through the assumption of local thermal equilibrium 
we find that at times tq ~ 0.6 fm cq — 15 — 20 GeV/fm^ at ^/snn = 200 GeV, well in excess of 
the critical energy for deconfinement (cc — 1 GeV/fm^). 

More detailed information about the properties of the quark-gluon plasma can be extracted 
via two particle correlation analysis. At present, a key question for perturbative QCD 
phenomenology is whether the medium induced gluon bremsstrahlung can significantly alter 
the di-hadron correlations measured at RHIC 123 ED- We naturally focus on the away-side 
|A(/?| > ^ case, where medium effects are the largest. A modihcation that does not change 
the A(/9-integrated cross section is vacuum- and medium-induced acoplanarity. The deviation 
of jets from being back-to-back in a plane perpendicular to the collision axis arises from the 
soft gluon radiation and transverse momentum diffusion in dense nuclear matter |22| . In the 
approximation of collinear fragmentation, the width of the away-side hadron-hadron correlation 
function can be related to the accumulated di-jet transverse momentum squared in the (/j-plane, 
sin Y^dpar = ^J^{k^)ip/p±d- Assuming a Gaussian form. 


( 7 ) 

a good description of |A(/3| > ^ correlations measured in elementary p-|-p collisions at 
Vs = 200 GeV dEI requires a large {k^ vac)‘P ~ ^ GeV^ for the di-jet pair with away- 
side scattered quark (and a 2.25 larger value for a scattered gluon). Additional broadening 
arises from the interactions of the jet in the QGP that ultimately lead to the reported 
energy loss, {k^y^V — 2 /Up)( 2 ;o)/\,g(- 2 o) In A, although only half is projected on the (/j-plane, 

{k^}^ = {kTvax)<P + likx hot)- 

Two competing mechanisms do, however, change the pt 2 dependence of the perturbative cross 
section, Eq. ®. First is the the parent jet “d” fractional energy loss e = AE^jEd^ which we here 
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for simplicity consider on average and evaluate by integrating Eq. ©• It leads to a rescaling of 
the hadronic fragmentation momentum fraction Z 2 ^ Z 2 !— e) |2l- If the energy loss is large, a 
second mechanism is invoked as a consequence. Hadronic fragments of the radiative gluons will 
increase the probability of finding low and moderate pt 2 particles associated with the interacting 
jet [ 23 . 

To calculate di-hadron correlations, we first map the jet structure of a hard 90*^- scattered 
parton on rapidity y ~ ry = — lntan(0/2) [6 being the angle relative to the beam axis) and 
azimuth (f>, 
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The approximately flat rapidity distribution of the away-side jet near 1/2 = 0 can be used to sum 
over all emission angles 9 G (^minj^max) C (0)"^) yielding 
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It is critical to note that projection on a plane coincident with the jet cone axis (the 0-plane in 
Eq. 0 is one example) efficiently masks the 0* —> 0 void in the angular distribution of medium- 
induced gluons reported in Eig. ^ Our conclusion is independent of the physical mechanism 
that depletes the parton (or particle) multiplicity in a cone around the jet axis. 

The end analytic result for the modification to Eq. © per average nucleon-nucleon collision 
in the heavy ion environment can be derived from the energy sum rule for all hadronic fragments 
from the jet, 
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Here, Zg = 

Numerical results, shown in the left panel of Fig.|21 correspond to triggering on a high px^ = 4 
- 6 GeV pion and measuring all associated 7 r'*“+ 7 r*^ + 7 r“. Depletion of hadrons from the quenched 
parent parton alone leads to a large suppression of the double inclusive cross section - a factor 
of 5 - 10 in central and a factor of 1.5 - 2 in peripheral reactions with weak pt 2 dependence. 

Hadronic feedback from the medium-induced gluon radiation, however, completely changes the 

( 2 ) 

nuclear modification factor It now shows a clear transition from a quenching of the away- 

side jet at high transverse momenta to enhancement at pt 2 + 2 GeV, a scale signihcantly larger 
than the one found in ESI STAR data in central and peripheral Au+Au collisions m is shown 
for comparison. 

Two-particle distributions in A+A reactions, calculated from Eqs. © and © for a pti = 
6 GeV trigger pion and two different pt 2 = Ij 4 GeV associated pions, are shown in the Fig. Q. 
Qualitatively, the medium-induced gluon component to the cross section controls the growth of 
the correlation width in central and semi-central nuclear collisions. Quantitatively, the effect 
should be even larger than the one estimated here, which is limited by the imposed 0 < 9* < ^ 
constraint. Experimental measurements of significantly enhanced widths for |A(/?| > ^ two- 
particle correlation in A+A collisions should thus point to copious hadron production from 
medium-induced large angle gluon emission. 

















2. Dynamical shadowing and energy loss in cold nnclear matter 

Single and double inclusive hadron production in reactions involving nuclei is also modified 
by dynamical shadowing effects. Suppression of the particle production rates arises from 
the coherent final state scattering of the struck small-x parton with several nucleons and the 
generation of dynamical parton mass Here, = 0.09 — 0.12 GeV^ for quarks 

is the scale of higher twist extracted from deeply inelastic scattering in heavy targets [24] . For 
electroweak vector meson exchange processes in DIS resumming the power corrections leads to 
an effective rescaling in the values of Bjorken-x. For p+A and A+A reactions the dynamical 
mass is amplified for the case of final state gluon scattering by the Casimir ratio CaICf = 2.25. 
It leads to significantly larger power correction effects in hadronic reactions. Effectively, 
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where the normalization relative to the proton is already included in Eq. (II 011 and Cd = 1 (2.25) 
for quarks (gluons) respectively. 

The left panel of Eig. 0 shows the nuclear modification RdA^PT) for hadron production in 
d+Au reactions m- We emphasize that the effect is large at small pT and relatively small 
at high pt- Note that at forward rapidity in Eq. m -i = p^. Data on forward y hadron 
production from BRAHMS is also shown [22]. The right panel of Eig. |3] presents the broadening 
and attenuation of the two particle back-to-back correlations. Such effects are found to be large 
at forward rapidity and at small pTi j PT 2 m- The coherent multiple scattering discussed here 
take place (in the c.m. of the collision) over a short interval t ~ RaIi Ra and for the case of 
A-|-A reactions precedes and complements the energy loss in the QGP which develops at time 
scales t Ra- 

In cold nuclear matter energy loss is anticipated to be significantly smaller when compared to 
the one in the quark-gluon plasma. It has been experimentally verified through measurements 
in d -|- Au reactions at RHIG at = 200 GeV where a small enhancement consistent 

with Cronin multiple scattering m has been observed m- It as been argued, however, that 
the radiative energy loss induced by the scattering of fast on-shell partons in nuclear matter 
evaluated by Bertsch and Gunion m can be implemented as a “Sudakov” suppression factor 
S{xf) ^ 1 — xp uni and can give a good description of the A“, a < 1 suppression in a large 
class of observed cross sections in hadronic reactions at forward xp- 


3. Conclusions 

In summary, we calculated the nuclear modification of inclusive moderate and high pp particle 
production and the back-to-back particle correlations in the framework of the perturbative QGD 
factorization approach, augmented by inelastic jet interactions in the quark-gluon plasma. Jet 
tomographic analysis of hadron attenuation points to initial energy density of 15 — 20 GeV/fm^ 
in ^sftn = 200 GeV central Au-|-Au collisions. At RHIG energies we found that the medium- 
induced gluon radiation determines the \A(p\ > ^ two-particle yields and the width of their 
correlation function to surprisingly high transverse momentum pp^ ~ 10 GeV. The predicted 
transition from back-to-back jet enhancement to back-to-back jet quenching is quantitatively 
compatible with the recent STAR data. In cold nuclear matter at RHIG coherent final state 
multiple scattering gives the dominant contribution to the attenuation of single and double 
inclusive hadron production. Energy loss, however, may also play an important role at 
forward xp- 
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Figure 4. Left panel: modification of the single inclusive hadron production from resummed 
nuclear enhanced power corrections in forward rapidity d+Au collisions at RHIC. Right panel: 
modification of the back-to-back di-hadron correlations versus rapidity, centrality and pT in 
y/SNN = 200 GeV d-l-Au reactions. 
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